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Abstract. We determine the critical coupling of the two-dimensional Brown- 
Ravenhall operator with Coulomb potential. Boundedness from below has 
essentially been proven by Bouzouina, whose work however contains a trivial 
error leading to a wrong constant exactly one half of the actual critical con- 
stant. Furthermore we show that the operator is in fact positive. Our proof of 
that is for the most part analogous to Tix's proof of the corresponding result 
for the three-dimensional operator. 



1. Introduction and partial wave analysis 
The two-dimensional free Dirac operator is defined in the Hilbert space L 2 (R 2 , C 2 ) 

by 

ch d ch d 2 
D = —ai- h —o- 2 - h mc o 3 , 

where h denotes Planck's constant, c the speed of light, m the electron mass and 
1 \ A) -A A 



'.I Oj' ° 2 = v z Oj' a3 -\0 1 

the Pauli matrices. D is a self-adjoint operator on H 1 (R 2 ,C 2 ), therefore we can 
define 

A+ := X(o,oo)(-D)- 

A + is obviously a projection and H+ := A + (L 2 (R 2 ,C 2 )) is called the positive 
spectral subspace of D. 

Finally the two-dimensional Brown-Ravenhall operator with Coulomb potential 
is defined in the Hilbert space H + by 

B = A + (VA) A+ , S>0 

By scaling this operator is unitarily equivalent to the one with m = h = c = 1 
up to the factor mc 2 and a change of 8 to 4-. In the following we will drop m, c 
and h. The positivity of the three-dimensional Brown-Ravenhall operator has been 
proven by Tix [B] and, by an alternative method, Burenkov and Evans 2;. Here 
we are interested in the two-dimensional result because of its importance in the 
description of graphene. 
Our main theorem is 

Theorem 1.1. Let 

Sr. : = 



r(i) 4 
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If 5 < 5 C , then 

B> (1-2.5), 

m particular B is positive. If 6 > 5 C , then B is unbounded from below. 
The approximate value of S c is 0.378. 

The hrst part of the proof has been carried out by Bouzouina \1\ , unfortunately 
with a small error that leads to a wrong value of 5 C . We proceed to sketch his 
development without proofs and correct the error where it was introduced. 

The free Dirac Operator corresponds in Fourier Space to the 2x2 matrix mul- 
tiplication operator given by 

D(p) =( 1 . Pl 

\Pl +V2 -1 

which has the pointwise eigenvalues ±E(p), where 



e(p) = V\pT+1- 

A + corresponds to the pointwise projection onto the eigenspace of E(p). Thus 
the Fourier transform of any element of % + can be written as u ■ £, where u is an 
element of I? (M 2 ) and £ denotes the normed eigenvector associated with E(p). So 
we can reduce the quadratic form belonging to B to one over L 2 (R 2 ): 

(i>,Bi>) = J E(p)\u(p)\ 2 dp-^f J u( P )^(p 7 )K(p,p') dpdp' =: (u,hu), 

(E(p) + l)(E(p') + l) + p& 



where 



K(p,p') = 



iV(p)iV(p')|p-p'| 



and 



N(p) := y/2E(p)(E(p) + 1). 
For u £ L 2 (M. 2 ) we can write 



u (p)=«(r e ") = ^£^=a fc (r)e iW , 
where the at are in L 2 ((0,oo)). Using this decomposition we get 



Theorem 1.2. For any u € L 2 (R 2 , y/l + \p\ 2 dp) holds 

(u,bu) = ^2(a k ,b k a k ), 

feez 

where 

(a k ,b k a k ) = e(r)\a k (r)\ 2 dr / / a k (r)a k (r')K k (r,r') drdr 

JO 7T Jo JO 

and for k > 

tf fe (r,r') = ft(ry)Q M/a Q (7 + 7)) + ^(r,r')Q fc+1/2 Q (7 + 7)) 
whereas for k < 

l/r r'\\ „ . _ fl / r r' 



K k (r,r') = ^(r,r')Q^ 1/2 Q (7 + 7)) + /3 2 (r ) r')Q_ fc _ 3 /2 Q ( 



r r 
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We have used 



(e(r) + l)(e(r') + l) 
Pi{r,r) = , w A , P2{r,r) 



n(r)n(r') ' n{r)n{r') , 



e(r) = \Jr 2 + 1, n(r) = v / 2e{r)(e(r) + 1). 

The Qi are Legendre functions of the second kind. 

Proof. This is Lemma 2.1 and Lemma 2.2 of pQ with the difference that the Kk (r, r') 
of [T] are incorrectly twice as big as here. In the proof of Lemma 2.1 Bouzouina 
claims that for any q g]0, 1[ and any 1,1' E Z 

2tt /.27T ue -We' f 2n cos (w) 

dOdO' = AttSu* I , H j d0, 



o 

when in fact 



y/l-qco&iO-P) ' Jo y/l-qcos(0) ' 



d8d8' = / e^-'')*' / 6 



o 7o v/l-gcos^-fl') 7o Vo v/l-gcos(0-0') 



as 

^ sm{W) 



VI — qcosO ' Jo VI ~ ?cos( 

d6> = 0. 



o VI - 9 cos ( 

□ 

The Legendre functions Qk-i/2i k <EN occurring in Theorem 1 1 . 2 1 are positive on 
]l,oo[ and (Qk-i/2)keN is a decreasing sequence (see [TJ Lemma 2.2]). From that 
it is easy to verify that K$ > Kk for all k G Z. Therefore we have 

Corollary 1.3. We have 

mf {(u,bu)\u £ L 2 (r 2 , VT+ |p| 2 dp) , ||u|| a = l} 
= inf {</,6 /)|/ € I 2 ((0,oo), vT+^dr) , ||/|| 2 = l} 
Proo/. This is Corollary 2.4 of [1]. □ 

2. POSITIVITY OF THE OPERATOR 

In view of Corollary [T73] it suffices to prove the equivalent of Theorem 1 1.1 1 for the 
form &o- 

Theorem 2.1. If 5 < S c , then 

(2.1) b > (1-2(5). 

Proof. The first part is an adaptation of the three-dimensional case in jB] . Therefore 
we will be brief. It suffices to prove (I2.1[) for S — S c . Then the general case is proven 
in the following way: 

(f,b f)=(l-£j J\(p)\f(p)\ 2 dp 

jr / POO £ pOO pQO \ 

e(p)\f(p)\ 2 dp- -^ / W)f(p)K (p,p')dpdp')> 



6, 



o Jo 



S. WALTER 



l~j-)\\f\\ 2 + j-(l-25 c )\\f\\ 2 = (l-'25)\\f\\ 2 
Finally we need only consider positive /. By a simple calculation one gets 



and 



2y e(p)Y e(p') 

By the usual application of the Cauchy-Schwarz inequality introducing positive trial 
functions ho and hi (for details see [5]) we have 



(/,6o/>> inf £(p) 

p€[0 : oo) 



£W;=e(p) -i,ff 1 + JL^ + A_ MWrt 



where 



2?r VV e b)/ fro(p) V e (p)/ fti(p) 

and 



I k {p) := / h k (j/)Q k - 1/2 [\(^ + -))dp'. 



, 2 VP' P 

It remains to choose ho and hi such that this infimum is 1 — 25 c . Note that choosing 
ho(p) — hi(p) = | as Evans et al. did in [3] would result in the infimum being 
finite, i.e. boundedness of the operator from below, but negative. Similarly to Tix 
we look for k = 0, 1 at the functions gk(re %e ) := r a e~ r e ike with a > —2, so that the 
gk are in L 1 (R 2 ). Let the Fourier transform on M 2 be defined as usual by 



/(P) := ^ / e- 4 P' x /(x) dx. 



It is a well known theorem [5J Theorem IV. 1.6] that for / £ L 1 (K 2 ) of the form 
f(re l6 ) = fo(r)e lke the Fourier transform is again of this form, more precisely: 



f(pe^) = e"*(-i)* / J k (pr)f (r)r dr, 
Jo 

where the Jk are Bessel functions of the first order. With [H 6.621.1] it follows that 

(2.2) g k (pe^) = e ik *(-i) k (p 2 + l)" 2 * 2 ^ + a + 2)P~ k 1 ((p 2 + 1)"*) . 

Let now a = —i. From the integral representation (4J 8.714.1] 

/2~ shr"(0) cos((i/+ eft , 1 

V7rr(i-^)J (cost - cos0) M+ 3 2 

and the fact that the gamma function is positive for positive argument it follows 
that the radial part of g^, i.e. 

Hp) := (p 2 + i)-t r (* + 1] p-f /2 (V + 1)* 
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is positive. Now we calculate the Fourier transform of Aofe in two different ways. 

\x\ 

On the one hand 

(n-)" ( ^-s(n-*)<^)-sjL^*'- 



_ e »**(_i)* / / p ' <ty'd p ' 
Jo Jo \J p 2 + p' 2 — 2pp' coscj)' 



dp', 



1 ( V , 

P , 

where the last identity follows from Lemma 2.2 of [1 . On the other hand, setting 
a = — | instead of a = — i in (|2.2[) , we get 



• <pe») = e«*(-i) V + l)"*r (* + i) (V + 1)"*) 



Now let 



hk( P ) := VP(^ + l)-*r ^ + - J Pf/* + 1)- 
which is positive. Then 



hU>) = K V Tir~ 1) ^(fc+^P_-* /2 (l/r f u 



Therefore 



and 



S{p) = e(p) - * c | (e(p) + i f-^ffi + (e(p) - 1) • | ■ -^TT 

inf £(p) = inf /(x), 
pG[0,oo) xe(o,i] 



where 



x 



Pi/2^) ' \x 7 3 p-i(x) 



Now we proceed to show that mf xe ro t -n /(x) = 1 — 2S C . According to 4, 8.704] 
we have 



r(l -/i) ^l-at 
and therefore 



X \Vx J P( 2' 2' S^) \ X ' ^ P ( 2'2'^' 1 2 X / 



By definition 



F{a, b; c; x) = > ,\ /■ x 



for |x| < 1, where 

(a)k := a{a + 1) • • • (a + n - 1), (a) := 1. 
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So we have for < x < 1 



F 



-v, v + 1; 1 — /x; 



1-x 



E 

fc=0 



{-v)k(v + l)fc 
(1 - M )fcfc!2 fe 



(l-z) fe . 



We see at once, that /(l) = 1 — 25 c . Examining the series above for v — ±.\ and 
[i = 0,-1 we notice that for v = — i all terms are positive, whereas for v = h all 
terms except the first one are negative. Furthermore we can estimate 



E 



(-t/)fc(t/+i) fc 

(l-ji) fc fc!2* 



fc=3 

So we have 



< 



(-v) 3 (v + l) 3 



(1 - /i) 3 3!23 



(l-a;) 3 -2 < 



(-i/)a(i/ + l)a 



(1-M) 2 2!2 2 



(1-x) 2 



/(s) - (i - 2<y > i - <y c ( i + 1 

x \x 



l + \{l-x) + £- % {l-x) 
l-Ul-x)-^{l-x) 



(1-2<5 C ). 



3 1-^(1 -x)- Ml'*) 2 
After some calculation this equals 

((49152-114688(5 c )(l-a;) + (-27648+48128 ( 5 c )(l-a;) 2 + (-4224+16128 ( 5 c )(l-x) 3 + 
+(1800 - 3504(5 C )(1 - a;) 4 + (225 - 5405 c )(l - xf)/ 
(3a;(128 - 48(1 - x) - 15(1 - x) 2 )(128 - 24(1 - x) - 5(1 - x) 2 )). 

The denominator is obviously positive and so are —4224 + 16128<5 C , 1800 — 3504<5 C 
and 225 — 540<5 C . Thus the term above is positive, if 

(49152 - 114688£ c ) + (-27648 + 48128<S C )(1 -x)>0 

and this is the case if x > 0.4. Therefore f(x) > (1 — 2S C ) for x > 0.4. 

To show that f(x) > (1 — 2S C ) for x < 0.4 we need the series expansion at 
[S 8.775]: 



2» cos {\{y + n)n) T (^±|±1) (j _ ^ g (f±f±i) fc 



+ 



2^+ 1 sin(l( t / + ^)7r)r(^ + l 



i 2 ; ta; 2fc + 



■x(l-x 2 f*Yl — 2 ' 



fe=0 



:§)**» 



for < a; < 1. So we have 

P-l/2 
1 



k y 



..2k 



„2fc+l 



fc) fc!' 



r(| 



2tt 



E 

fc=0 



-1/2 



(x) 



r(| + fc)r(-j + fe) afc 
r(i + fc)fc! 



fc) fc! y 

r(i + fc)r(| + fc) 
r(| + fc)fc! 



,2fc+l 



47T 



fe=0 



r(i + fc)fc!" r(f + fc)fcf 



„2fc 



„2fc+l 



vi ^i /2 W- 47r 2^ r(i + fc)fc! 



.2 A' 



fc=0 



r(| + fc)r(-i + fc) ^ 2fc+1 



r (§ + *)*! 



POSITIVITY OF THE TWO-DIMENSIONAL BROWN-RAVENHALL OPERATOR 

To arrive at an estimate we first show that the function mapping x to „F^ X \ \ 
strictly decreasing on (0, oo). For x > 

d T(x) _ r(x)r (x + |) - r(x)T' (x + |) _ r(a:)V(a) - r(a:)V (a + |) 



is 



da; r (a; 
where 



Tlx 



r(x 



tp{x) := 



T(x)- 



According to 01 8.361] 



dt 



and therefore ip is strictly increasing and 



r(aQ 



is strictly decreasing on (0, oo). 



If we denote the nth coefficient of the respective power series by a n , we get in each 
of the four cases for k > 1 



0>2k 



and 



02fc+l 



> 



r(| + fc)r(i + fc)r(| + fc) _ fc + 



»2fc+2 



r(| + fe)r(| + fe)r(i + fe) k + \ 
r(| + fc)r(f + (fc + i)! fc + i 



> i 



r(| + fc)r(| + fc) fc! 



> i. 



So starting from a2, the coefficients alternate and their absolute values decrease 
strictly. Therefore 



P-1/2O) < 
Pi/i(x) > 



y/l-x*PZi /2 {x) < 




y/l - x 2 P-/ 2 (x) > 

One can easily check that each of the right sides is positive for x < 0.4. Plugging 
this in and cancelling we get 

1 /rfi) 4 ■ 1 

m - { l-26 c) 



8tt 2 

r(i) 4 




8tT 2 - ^(^J + TrV 

i2^-|r(i) 4 



r(i) 4 -4A + ir(i) 4 ^ 
871-2 + |r(i) 4 x-37r 2 x 2 
r(i) 4 Stt 2 



r(!) 



8tt 2 



r(i) 4 



The right side equals 



(r ( i) 8 + 64tt 4 ) (8r (i) 4 - 96^ 2 a; - 3r (±) 4 x 2 ) (i6tt 2 + T (±) 4 x - 6ir 2 x 2 
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where 

b o = 8 T ( =- } - 256 T ( - A I tt' z + 3072 r ( - A ) tt 4 - 98304tt s 



1\ 16 „ /n 12 , /l^ 8 



b 1 = -8v(^j + 3072 rQ^) tt 4 - 40960 r tt 6 + 98304tt 8 

&2 = -3rQ^j + 192rQ^) 7r 2 -2944rQ^ 7r 4 + 4096 rQ^ 7r 6 + 368647r 8 
b 3 = 3rf^\ -24T^~) 7r 2 -128rQ^ tt 4 + 10752 T (jj tt 6 - 36864tt 8 

b 4 = -12T (j\ 7r 2 + 288 rQ^ 7r 4 - 1536 rQ^ tt 6 

One easily verifies that in the fraction above the denominator is positive as well as 
bo and 03, while b\, 02 and 04 are negative, so 

4 

h nx n >b -h- 0.4 - b 2 ■ OA 2 - b 4 ■ 0.4 4 > 0. 



n=0 



3. Unboundedness above the critical constant 
The following is inspired by [3]. 
Lemma 3.1. Let v = ±1 and 



□ 



There exists C € R suc/i i/iat /or all \ < a < b 

fjf^«"(KM))**"-*G) +c 

Proof. According to [7, p. 334] for i > 1 and v > -1 we have the integral represen- 
tation 

.t-VF^i M t -s/W=I u-1 

Quit) = / - ; < (t- V* 2 -l) / ; ' dx. 

W 7o Vx 2 - 2ta + 1 " V Vo Vx 2 - 2ta + 1 



If < p < 1 and i = | + ij , then t — \Jt 2 — 1 = p; we conclude that 
According to [4j 8.832.3] the derivative of the Legendre functions is given on (1, oo) 

by 

d . (i/ + 1)(Q 1/+ i(i) -a; Q v (x)) 

= ■ 

Let a € R and < p < 1: 
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) 








O-K 1 -*) 


h 


j> 2 




i 

2 



2(- + 1) (i (p + 1) Q, (*(* + ?))" Qh-i (I (f + ?))) • 



2(^ + l)(l(p+i)Q,(l(p+i))- P Q, (1(^+1^ 

= (i/ + l + a)Q„Q 
Therefore the function 



„a-l 



< 

and 

r 6 



2 V P/ 

is increasing on (0, 1), if a > — v — 1, which is satisfied by f = ±\ and a = — \. 
Thus 

jf;0*(^^))*'*-jr;j( S y*(i(^y))*'* 

d 

r;fMKH))*^i)i s MK>'v))^ 

£ 2 1 ? *■ (K? + ;)) * - 2 f *■ G (' + ;)) * - 

From this the assertion of the lemma follows with C — —11 v . □ 
Theorem 3.2. If 8 > 5 C , b a is unbounded from below. 
Proof. Let 1 < a < b and define for p > 

/(p) := X(o,6)(p) ' ~ v 

We have for p,p' e (a, 6) 

/Mp,p') > \ 

and 



2 V e(p)V e(p') " 2 \ e{p)J\ e{p') 
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1 1 



>rU j 1 > 

P P 



1 1 



From 1 we know that 



I-i 



/2 



47T 



and I 1 / 2 



16tt 3 



The facts collected above and Lemma I37L1 arc used in the following estimate. From 
now on we replace some terms that do not depend on a or b by the word 'const': 



~ / f ~ ~J k o(P,p') dpdp' 

I" J a J a P P 




const 



const 



const. 



Furthermore 



Thus 



/ e(p)~2 rfp < / — dp < In I ~ 



(f,b f)<c 1-- 1-- In - 



const. 



Since j- > 1, we can choose a large enough so that 



< 0. 



1-^11-- 

a 



If we keep a fixed and let b go to infinity, (/, bof) tends to negative infinity, while 
1 1/| | 2 = (I - I) < 1. This concludes the proof. □ 

The proof of Theorem If .H is now complete. 

Acknowledgement. The author thanks H. Siedcntop for giving valuable advice. 
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